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Abstract 

The peak of the de-absorbed energy distribution of the TeV emitting blazars, all 
of the BL Lacertae (BL) class, can reach values up to ~ 10 TeV. In the context 
of synchrotron-self Compton (SSC) models of relativistic uniformly moving blobs 
of plasma, such high energy peak emission can be reproduced only by assuming 
Doppler factors of 5 ~ 50. However, such high values strongly disagree with the 
unification of FR I radio galaxies and BLs. Additionally, the recent detections of 
slow, possibly sub-luminal velocities in the sub-pc scale jets of the TeV BLs MKN 
421 and MKN 501 suggest that the jets in these sources decelerate very early to 
mildly relativistic velocities (T ~ a few). In this work we examine the possibility that 
the relativistic flows in the TeV BLs are longitudinally decelerating. In this case, 
modest Lorentz factors (T ~ 15), decelerating down to values compatible with the 
recent radio interferometric observations, can reproduce the ~ few TeV peak energy 
of these sources. Furthermore, such decelerating flows are shown to reproduce the 
observed broadband BL - FR I luminosity ratios. 
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1 Introduction 

The blazar TeV emission reaches the Earth partially absorbed by the Diffuse 
InfraRed Background (DIRB) radiation which pair-produces with the TeV 
photons (|l|). The de-absorbed spectra depend on the source redshifts and 
the still elusive energy distribution of the DIRB. However, for all expected 
DIRB forms, both the intrinsic peak energy E p and peak luminosity L p of the 
TeV spectral component are higher than those observed. Even for the nearby 
(z = 0.031) MKN 421, E p can increase by a factor of ~ 10 after de-absorption 
to ~ a few TeV (2). The de-absorbed spectrum of H1426+428 at z=0.129 is 
even more extreme, characterized by E p ~ 10 TeV (J3|). 
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The BL synchrotron spectra exhibit a break at energies e' b ~ 10~ 4 — 10~ 6 
(primes denote energies on the flow rest frame, all energies normalized to m e c 2 ) 
such that the greatest fraction of the synchrotron luminosity is at energies 
e' > e' b . For this reason and because of the K-N decrease in the Compton 
scattering cross section, electrons with 7 £ l/e' b will inverse Compton scatter 
a decreasing fraction of the available photons; as a result, the maximum IC 
luminosity will occur at energies e' p ~ or e p e b ~ 1. The values of e p , e' b 
observed at the lab frame are e p = 5e' p and e?, = 5e' b , yielding the following 
relation between 5 and the observed energies e p , €b 

5 £ (e 6 e^) 1 / 2 = 40 (i/ 6>16 ^,ioTev) 1/2 , (1) 



where v bj i$ is the observed synchrotron break frequency in units of 10 16 Hz, 
E p , iotcV is the energy of the de-absorbed TeV peak in units of 10 TeV and 
5 = l/r(l — (3 cos 9), with j3 the dimensionless flow speed and 6 its angle to 
the observer's line of sight. These values of <5(~ Y) are in strong conflict (0) 
with the unification of BLs and FR I radio galaxies (0), which require Lorentz 
factors r rsj 3 — 7. Additional constrains for the inner jet flow of the TeV 
blazars come from the small, possibly sub-luminal apparent velocities observed 
interferometrically in the sub-pc scale jet of MKN 421 and 501 (0; S 0), 
suggestive of a decelerating flow in the inner jet of TeV blazars (j3). 



2 Decelerating relativistic flows: application in TeV blazars 



Motivated by the above analysis and the unification of the multiwavelength 
spectra of the hotspots of FR II radio galaxies and quasars in terms of de- 
celerating relativistic flows (10), we propose that the same considerations are 
applicable in resolving the conflict of the high <5's with TeV blazar unification. 
So, we assume that a power law electron distribution is injected at the base of 
a relativistic flow that decelerates while at the same time cooling radiatively. 
The highest frequencies of its synchrotron component originate preferentially 
at its fast base where the electrons are more energetic and its Lorentz factor 
largest. As both the flow velocity and electron energy drop with radius, the lo- 
cally emitted synchrotron spectrum shifts to lower energies while its beaming 
pattern becomes wider. The observed synchrotron spectrum is the convolution 
of the comoving emission from each radius weighted by the beaming amplifi- 
cation at each radius. At small angles the observed spectrum is harder than 
that observed at larger angles. This is the result of the progressively smaller Y 
of the flow with distance in combination with the concomitant lower electron 
energies due to cooling. 

The inverse Compton emission of such a flow behaves in a more involved 
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Fig. 1. Left panel: the Synchrotron and IC emission from a decelerating relativistic 
flow r = 15tor = 4at(9 = 3° (solid line) and 9 = 6° (broken line) (see text for 
details). The broken vertical line corresponds to the maximum IC energy ^ max 5 for 
6 = 3°. The shaded areas correspond approximately to the energy range of X-ray 
and TeV telescopes and the two vertical solid lines energies to 1 and 10 TeV. Right 
panel: The SED of a decelerating flow for a range of observing angles. The physical 
parameters are the same as before, but with ^y-max = 2 x 10 5 . The shaded bars reflect 
the average luminosity range in radio, optical and X-rays, between the samples of 
BLs and FR I radio galaxies (|lll). 

way: Electrons will upscatter the locally produced synchrotron seed photons, 
giving rise to a local SSC emission with 5— dependence similar to to that of 
synchrotron. However, the electrons of a given radius scatter not only the 
locally produced synchrotron photons, but also those produced downstream 
in the flow. The energy density of the latter, will appear Doppler boosted in 



the fast (upstream) part of the flow by ~ T^ e/ (|12l). where T re[ is the relative 
Lorentz factor between the fast and slow part of the flow. With their maximum 
energy being lower (because of cooling) and their energy density amplified they 
can now contribute to the IC emission at energies higher than expected on 
the basis of uniform velocity models (see section 1). Also the 5— dependence 
of this upstream Compton (UC) emission will be different from that of SSC 
and more akin to that of external Compton (EC) (fl2Tb 

In the left panel of fig. 1 we plot the SED for a flow decelerating from T = 15 
to T = 4 for two angles 9 = 3° and 9 = 6° over a distance Z = 2 x 10 16 cm. The 
radius of the cylindrical flow is set to R = Z = 2 x 10 16 cm, while the electron 
distribution at the base of the flow is 71(7) oc 7 -2 , 7 < 3 x 10 7 and the magnetic 
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field B = 0.1 G, half its equipartition value. At 9 = 3° this model achieves 
a peak energy for the high energy component at ~ 10 TeV, using a modest 
initial Lorentz factor of T = 15. Note the strong dependence of synchrotron 
peak energy on 9 discussed above. Finally, note that, in contrast to uniform 
velocity homogeneous SSC models, the Compton component is more sensitive 
to orientation than synchrotron, an indication that UC scattering dominates 
the ~ TeV observed luminosity. 

We now turn our attention to the problem of the unification of BLs with 
FR I sources. A comparison of a sample of FR I nuclei to BLs of similar 
extended radio power shows that FR I's are overluminous by a factor of 10— 10 4 



compared to their expected luminosity (|5tllll). under the assumption that BLs 
are characterized by Lorentz factors T ~ 15, and that FR I's are seen under an 
average angle of 60°. In particular the average BL to FR I nucleus luminosity 
ratio at radio, optical and X-ray bands was found to be: \og(L BL / L FRI ) R pa 
2.4, \og{L BL /L FRI ) opt pa 3.9, \og{L BL / L FRI ) X ~ 3.5. In the right panel of 
fig. 1 we plot as vertical bars the luminosity separation of BLs and FR Is 
according to (jslllll). We also plot the SED of a decelerating flow with physical 
parameters similar to the one described above but with 7 max = 2 x 10 5 to 
produce an SED synchrotron peak similar to those of the intermediate BLs 



that correspond in extended radio power to the FR Is of (J5J; lllf ). As can be 
seen, the range of the model SEDs between 9 = 60° (FR I) and 9 = 1/T (BL) 
reproduce well the observed range in luminosities, while a uniform velocity 
model of T = 15 (as demanded by fits of the 7— ray spectra) would produce at 
9 = 60° an SED many orders of magnitude smaller than shown in the figure. 
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